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°^ ' Abstract 
m 

On . 

A higher-twist nuclear enhancement of R = o~l/o~t, as might be expected 

m 

to arise due to fermi motion and whose magnitude is within the error-bars of 
recent experiment, is shown to lead to a monotonic decrease in the ratio of 
nuclear vs. nucleon cross-sections at small Bjorken x for increasing Q 2 . This 

• 

effect at small x, comparable in magnitude to those reported for shadowing, 

is driven mainly by kinematic factors and essentially vanishes for x > .1. Its 

unusual Q 2 dependence rather complicates the unravelling from present data 

in the shadowing region the corresponding dependence in Q 2 of the nuclear 

structure functions, F^(x,Q 2 ). 
PACS numbers: 25.30.Mr, 13.60.Hb, 12.38.Qk 
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Shadowing, the observed depletion for heavy nucleii in the nuclear vs. nucleon cross- 
section at low Bjorken x in real or virtual photon scattering, is a subject of both continuing 
experimental and theoretical interest. The experimental situation can roughly be seperated 
into two categories: those at low (or zero) photon virtuality, Q 2 , and those at moderate 
to large ||, H values of Q 2 . It is the data at large Q 2 that is of concern for this study, for 
which Ref. || is of particular interest as it is there that the double-differential cross-section 
d 2 a/dxdQ 2 is presented for a range of Q 2 for a number of x values in the shadowing region. 
What is to be here addressed is our ability to infer from this data the correct Q 2 dependence 
of the structure functions F 2 A (x,Q 2 ) in the shadowing region. This is important not only 
for hopefully distinguishing and potentially selecting different theoretical models [|J, || for 
shadowing at large Q 2 , but because it is also crucial for our understanding the results of 
other experiments such as Drell-Yan |J or J/ip production J7|, in which other mechanisms 
||, might be playing an important role. It should be emphasized at the outset that this 
is not an attempt to explain shadowing, for which a diverse set of models presently exist 
and which can be roughly seperated into those using quark and gluon H and those utilizing 
hadronic degrees of freedom. Indeed, the effects here discussed all but vanishes for the 
smallest Q 2 values of Ref. Instead it is the point of this paper to demonstate how a 
very simple, higher-twist nuclear effect (i.e. which vanishes with increasing Q 2 ) leads to a 
rather unusual Q 2 dependence in the double-differential cross-section, and one which would, 
if verified, greatly complicate our ability to infer F^(x, Q 2 ) from existing data. 

To order a 2 , the double differential cross-section per nucleon for lepton scattering from 
an unpolarized nuclear target may be written as: 

d 2 a _ A7ra 2 e F 2 (x,Q 2 ) 
dxdQ 2 Q 4 x 

where M = .94GeV is the mass of a (free) nucleon, E is the incident lepton's en- 
ergy, F 2 (x,Q 2 ) is the nucleon's structure function (in a nuclear medium), and R = 
ol/&t is the ratio of longitudinal to transverse cross-sections, which in terms of the 
usual functions F ± (x, Q 2 ) and F 2 (x,Q 2 ) is given by R(x,Q 2 ) = F L (x,Q 2 )/2xF 1 (x,Q 2 ), 
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Q 2 Q 2 Q 4 f l+4M 2 x 2 /Q 2 \ 

2MEx AE 2 + 8M 2 E 2 x 2 \ 1 + R(x, Q 2 ) ) 



where F L (x,Q 2 ) = (l + ^t~) F 2 (x,Q 2 ) — 2xFi(x,Q 2 ) is the so-called longitudinal struc- 
ture function. It is well known that upto target mass corrections, R vanishes in the parton 
model for the case of spin 1/2 partons due to helicity conservation. In perturbative QCD, 
R acquires a non-zero value from a s (Q 2 ) corrections to the structure functions. In addition 



to these corrections, R is also non-zero due to an explicit higher-twist correction [ID|, [IT 
to F L (x,Q 2 ): 
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4A 

FL = \x,Q 2 ) = —T 1 (x), (2) 

in which T\(x) is a new nucleon matrix element. It involves transverse components of the 
covariant derivative acting upon quark field operators. If one ignores the fact that this is 
a covariant derivative (required for gauge-invariance) , this new matrix element corresponds 
naturally to what one associates with an intrinsic transverse momentum of the quarks in 
the parton model. That is, one obtains in the limit that the gluon coupling in the covariant 



derivative g — ► 0: [fLlJ 



, 4A 2 T!(:r) 4 r f2 ,. , 2 . (kV) 

lim ^ = — J d 2 k ± k 2 J(x, k\) = 4^, (3) 



g-+o Q 2 Q 2 J - ^ V ' ^ Q 

where in the parton model f(x, kjj is related to F 2 (x) by 



Fl= 2 



(x)=xjd 2 k ± f(x,kl). (4) 



Recent theoretical fits [|T^] to R for the free nucleon indicate that such a higher-twist element 
is indeed required to fit the data. 

If we now turn to nuclear targets, it is clear form Eq. ([I]) that a nuclear dependence 
to R could distort at any x or Q 2 the extraction of F 2 A {x,Q 2 ). Indeed, such a conjectured 
dependence |16j was in fact proposed to account for the initial discrepancies between the 
EMC [|I7]] and SLAC [18j data for x ~ .1. Besides the ultimate experimental resolution of 
this discrepancy, the shift in nucleii used by AR w .15, (taken to be essentially Q 2 



independent, as would be neccessary to fit the wide range of Q 2 available at these large 



values of Bjorken x) does not seem to be supported by later experiments |1J| at SLAC 



devoted to measuring AR. We're now though interested in AR at small x, for which Ref. 



j~5| is the sole source of experimental information. This data (as was that from |L4]]) is 
consistent with a value for AR = 0, but with such large error bars that the best one should 
probably safely conclude is that AR < .1. However it should be emphasized that because of 
the experimental difficulty in measuring this quantity, a systematic double-binning in both 
x and Q 2 for AR(x, Q 2 ) in either [|l4j or is lacking. 



A simple reason to expect a nuclear dependence for R is fermi motion. Returning to the 
heuristics of the partonic model, motion of the nucleons in the nucleus increases the average 
transverse momentum of the quarks inside a nucleon. By Eq. (|3|), such an origin for AR is 
higher-twist, and hence vanishes as Q 2 — > oo. Despite this Q 2 dependence, it will be seen 
that the effect on the measured cross-section d 2 a / dxdQ 2 at given fixed x actually grows with 
Q 2 . This growth is a direct consequence of the kinematic factors entering Eq. ([]]), as is also 
the fact that for a given fixed Q 2 , the effect grows for decreasing x. 



The world's data for R for a free nucleon has been collected and analysed by Ref. [13 
for which an empirical best fit has been obtained. It is given by the formula: 

ryfit = h ( Q 2 0.125 2 \ h b 3 

ln(Q 2 /A 2 ) V Q 2 + 1 0.125 2 + xV Q 2 Q 4 + 0.3 2 { ) 

where A = .2GeV, = {.0635, .5747, —.3534} and all momenta are in GeV. Note that 
there is a misprint in Jl3| for the value of the first coefficent, b\\ the decimal point needs to 
be shifted one place to the left .Q Whereas the expression multiplying the inverse-log term 
has been chosen to fit the results from pQCD for R at large Q 2 , one should not interpret too 
seriously the various coeflicents {6j} in terms of a twist expansion, as target mass corrections 
must first be disentangled. Such a seperation was undertaken in Ref. ||12j| . For our purposes, 
Eq. (H) is used to merely generate the best information we presently have on R for a free 



1 This misprint might account for the apparent discrepancy between certain statements of Ref. 
U and the results reported here as to the relative importance that a nuclear dependence to R 
generates in the differential cross-section. 



nucleoli. One caveat is that no data actually exists for R for x < .1. This fact will be 
returned to later. 

A higher-twist nuclear enhancement to R^ lt is conjectured. To R^ lt is added a term AR, 

R A = R fit + AR ^ (g) 

modelled to fall as 1/Q 2 ala Eq. @ with k± = .3GeV ~ kf ermi taken. No x dependence 
is given this enhancement although such a dependence is likely, especially in the small x 
region where the heuristic interpretation of the higher twist matrix element T\(x) is likely 
to break down ||11|| . In this regard, the simple connection with fermi motion is given merely 



as a plausibility arguement and the skeptical reader could interpret the following as simply 
an exercise exploring the effects that introducing such a higher-twist nuclear dependence 
for AR makes on the differential cross-section. 

Fig. [I] plots the resulting ratio of nuclear to nucleon cross-sections, Eq. |], such a AR 
generates for various values of Bjorken x ({xi} = {.0125, .025, .05, .1}) assuming no change 
in the nucleon to nuclear structure functions, i.e. F^ = F^. The kinematics is that of Ref. 
H (E = 200GeV) for which the curves cutoff for each at the kinematic upper limits of 
the experiment. From top to bottom the curves correspond to decreasing X{. Note the trend 
with Q 2 for each curve. While the plots in Fig. [I] are all consistently less than the shadowing 
actually seen, especially at small Q 2 (hence the initial statements that this is not an expla- 
nation of shadowing), for the larger Q 2 values the suppressions are large enough to play an 
observable role (e.g. the shadowing for Calcium at x — .025 is roughly 5% j3]). Observe 
most importantly that the effect at fixed Q 2 grows for decreasing x. This is dramatically 
demonstrated by the plot for x = .0125, where though it is admittedly a point for which 
both R* %t and AR might be most seriously doubted. It has been included nevertheless for 
emphasis. In Fig. the solid line is the x independent nuclear-enhancement, AR, added to 
Rf lt (x, Q 2 ) of Eq. ([5]). For reference is also plotted R^ lt (x, Q 2 ) for two illustrative values of 
Bjorken x. The dashed line is R^ lt (x, Q 2 ) for x = .1 and the dotted line is the (interpolated) 
value for x = .025. 



The results in Fig. [l| are significantly more sensitive to the assumed size of AR in the 
shadowing region than they are to R itself. For example, arbitrarily multiplying W u by 
a factor of two makes a change of only about .5% and 1% in the ratio of cross-sections 
for x = .025 and x = .0125 respectively. Since R for the nucleon is unknown in the 
shadowing region, this insensitivity is relevant. On the other hand, doubling (or halving) 
AR effectively doubles (or halves) the effect at large Q 2 for these same x values. As for the 
actual magnitude of AR being introduced, for Q 2 > 5GeV 2 , AR < .065 and would appear 
to be small enough to fit within the uncertainties (at any Bjorken x) of experiment fl4|], |15 



Indeed, for Q 2 = 9GeV 2 , where the effect in Fig. [T] for x = .025 is largest, AR < .04, well 
within experimental error bars. For smaller values of Q 2 however, the x independent value 



of AR is in apparent contradiction with the measured values [|14j] at x > .1. While this fact 
only further emphasizes the result contained in Fig. [I] that a nuclear dependence for R yields 
negligible effects on the measured cross-section outside of the shadowing region, theoretically 
one notes that for the smaller Q 2 values, the magnitude of the AR being introduced appears 
less and less as a perturbation to R(x,Q 2 ). For example, at Q 2 = 2GeV 2 , the used value 
for AR = .18 would be over a 60% correction to R(x, Q 2 ) at x = .1. In such circumstances 
corrections of order Q~ 4 cannot be justifiablely ignored, and indeed a twist expansion may 
not even converge fl9j . For the future experiments at CEBAF (E = 6GeV) that will measure 



AR(x, Q 2 ), this last consideration suggests that only at the very highest Q 2 there accessible 
might a twist-expansion in fact appear. For of course the main focus of this study, the 
shadowing region, a machine with much higher energy (and presumably an electron beam 
for higher intensities in order to improve upon the CERN data, as in e.g the proposed 50GeV 
SLAC upgrade) is neccessary for probing AR(x,Q 2 ). 

In conclusion, the Q 2 dependence of Fig. [I] is rather startling. In a counter-intuitive 
fashion, an intrinsically higher-twist effect for AR is of growing importance with increasing 
Q 2 in the observed rate. Examining Eq. (|IJ) this apparently surprising dependence with 
Q 2 occurs because near the ends of phase-space, 1/xE enters as an important new scale 
where the parameter y = Q 2 jlMxE (the fractional energy loss of the electron) is reaching 

6 



its kinematic upper bound: y < 1. In support of this effect, a close scrutiny of the data of 
Ref. H does indicate an apparent systematic increase with Q 2 of the shadowing observed 
for most of the values of Bjorken x < .1. While a nuclear enhancement of R would naturally 
account for such a trend in the data, such an enhancement would also, depending on its 
size for any particular x in the shadowing region, act to mask the true Q 2 dependence of 
F^(x,Q 2 ). Without therefore better knowledge of AR(x,Q 2 ) in the shadowing region, the 
extrapolations of present data to either yet higher Q 2 or other processes must be considered 
problematic. 
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FIGURES 

FIG. 1. The ratio of nuclear to nucleon differential cross-sections generated by a 
AR oc 1/Q 2 , as a function of Q 2 for given x. From bottom to top, the curves correspond to 
{xi} = {.0125, .025, .05, .1}. 

FIG. 2. In solid line, the value of the nuclear enhancement Ai? being introduced. For reference, 
also the values of R^ lt {x = .1) in dashes, and R^ lt {x = .025) in dots. 
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